The effect of initial fluctuations on bottomonia suppression in relativistic heavy-ion 

collisions 
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Using the screened Cornell potential and the next-to-leading order perturbative QCD to deter- 
mine, respectively, the properties of bottomonia and their dissociation cross sections in a quark-gluon 
plasma, we study in a 2+1 ideal hydrodynamics the effect of initial fluctuations on bottomonia pro- 
duction in relativistic heavy-ion collisions. We find that while initial fluctuations hardly affect the 
yield of the IS ground state bottomonium, their effect on that of excited bottomonium states is not 
small. Compared to the case with smooth initial conditions, the survival probability of excited bot- 
tomonia is reduced at low transverse momentum and increased at high transverse momentum. The 
observed suppression of the excited bottomonia relative to the ground state bottomonium by the 
Compact Muon Solenoid (CMS) collaborations at an average transverse momentum can, however, 
be described at present with both smooth and fluctuating initial conditions. 

PACS numbers: 



I. INTRODUCTION 

Since the suggestion of J/ip suppression as a possible 
signature of the quark-gluon plasma (QGP) formed in 
relativistic heavy ion collisions [l| , extensive studies have 
been carried out both theoretically and experimentally 
not only on charmonia production but also on bottomo- 
nia production in these collisions 0-fHl- Compared to 
J/ip production, the ground state bottomonium is, how- 
ever, a more promising probe to the hot dense matter 
created in relativistic heavy ion collisions because of the 
relatively small contributions from its excited states and 
from regeneration in the QGP [TU, [22|. Recently, suppres- 
sions of bottomonia production in relativistic heavy-ion 
collisions relative to those expected from p+p collisions 
at same energies have been observed at both the Rel- 
ativistic Heavy-Ion Collider (RHIC) [H| and the Large 
Hadron Collider (LHC) [IE El- To understand these 
experimental results, we have studied in Ref.(l7j bot- 
tomonia production in these collisions by including their 
production from both initial hard collisions of nuclcons 
and the regeneration in the produced quark-gluon plasma 
based on a schematic hydrodynamics for the bulk colli- 
sion dynamics [l7j . We found that the contribution from 
regeneration was very small for bottomonia as expected 
and the modification of the thermal properties of bot- 
tomonia in hot dense matter was helpful in describing 
the experimental data. The above study has, however, 
neglected the temperature fluctuation in the produced 
hot dense matter as all thermal quantities are taken to 
be uniform in the schematic hydrodynamics model. In 
the present study, we improve our previous results by 
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using a 2+1 ideal hydrodynamics to include also the ef- 
fect of initial fluctuations on bottomonia production in 
relativistic heavy ion collisions. 

The paper is organized as follows: In Sec. [Hj we intro- 
duce the 2+1 ideal hydrodynamic model for relativistic 
heavy ion collisions with both smooth and fluctuating 
initial conditions. We then describe in Sec. Illll thc prop- 
erties of bottomonia in QGP based on the screened Cor- 
nell potential for determining their dissociation temper- 
atures and the next-to-leading order perturbative QCD 
(pQCD) for determining their decay widths. Results and 
summary are given in Sec. IIVI andfVl respectively. 



II. THE 2+1 IDEAL HYDRODYNAMICS FOR 
RELATIVISTIC HEAVY ION COLLISIONS 

The hydrodynamic equations are based on the conser- 
vations of energy-momentum and various charges. Be- 
cause chemical pot entials are negligible in heavy-ion col- 
lisions at LHC |18| . only the energy- momentum conser- 
vations are used in this study. Assuming the boost in- 
variance, energy-momentum conservations are then ex- 
pressed in the (r, x, y, 77) coordinate system, with r = 
Vt 2 - z 2 and r) = \ In f±§, as 0, [U 

d T (TT 00 ) + 8 x (tT 0x ) + d y (rT°y) = -p, 
d T (TT 0x ) + 8 x {tT xx ) + d y {rT xy ) = 0, 



d T {TT 0y ) + d x {rT xy ) + d y (TT yy ) 



0, 



(1) 



where andp are, respectively, the energy-momentum 
tensor and pressure. 

To solve Eq. ([T]) requires information on the initial con- 
ditions of a collision, particularly the initial entropy den- 
sity, and the equation of state of the produced matter. 
For the initial entropy density, it is taken as 
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In the above, n par t and n co ii are, respectively, the number 
densities of participants and binary collisions. In the case 
of smooth initial conditions as used in Ref. 0, they are 
given by 
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where To is the initial thermalization time, which is taken 
to be 1.05 fm/c as in our previous study @, [Hj]; A = 
B = 208 is the mass number of Pb; r and b are, respec- 
tively, the transverse position vector and impact param- 
eter; (Ti n = 64 mb is the nucleon-nucleon inelastic cross 
section for LHC energies (24|; Ta(b) = / dzp A (B)(fi z ) is 
the thickness function with p A (B) being the nucleon dis- 
tribution function in nucleus A{B) for which the Wood- 
Saxon model is used. 

In the case of fluctuating initial conditions, the posi- 
tions of colliding nucleons are determined according to 
PA(B) by the Monte Carlo method. If the transverse dis- 
tance between a nucleon from nucleus A and a nucleon 
from nucleus B is shorter than \J a-m/ir, the two nucleons 
are then considered as participants and a binary colli- 
sion takes place at their middle point. In this case, the 
number densities of participants and binary collisions are 
given, respectively, by 
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where rj and fj are transverse positions of participant 
i and binary collision j, respectively. Here we use the 
same smearing parameter a = 0.4 fm for both number 
densities [2o| . 

The parameters C and a in Eq. ([3]) are determined from 
fitting the centrality dependence of the charged-particle 
multiplicity [21] . Using the Cooper- Frye freeze-out for- 



mula and assuming that the multiplicity does not change 
after chemical freeze-out at temperature T = 160 MeV 
[IH, we obtain C = 24.8 and x = 0.15. We note that this 
value of C is slightly smaller than that in our previous 
study based on a schematic hydrodynamics [22j . 

For the equations of state, we use the quasiparticle 
model based on the lattice QCD data for the QGP and 
the resonance gas model for the hadron gas as in Refs. @, 
|26| . This model thus assumes the presence of a first-order 
phase transition and the critical temperature T c is 170 
MeV. 

We solve the hydrodynamic equations [Eq. ([T])] numer- 
ically by using the Harten-Lax-van Leer-Einfcldt (HLLE) 



algorithm [27H29j . The accuracy of the calculation can 
be monitored by the entropy conservation condition , 



dS t , 



di] 



dxdy tsj±, (5) 



where s is the entropy density and 7^ = (1 



-^)- 1/2 

with v x and v y being the transverse components of the 
flow velocity. It is found that the entropy is conserved 
within 2 % in the case of smooth initial conditions, while 
it increases to 5~6 % in the case of fluctuating initial 
conditions. In the latter case, we thus rescale the entropy 
at each time step accordingly to restore the conservation 
of entropy 
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FIG. 1: (Color online) Initial temperature distributions in the 
transverse plane at to for smooth (upper panel) and fluctu- 
ating (lower panel) initial conditions in Pb+Pb collisions at 
^/sjviv = 2.76 GeV and impact parameter b = 2.1 fm. 

In Fig. [TJ we compare the temperature distributions 
in the transverse plane at the thermalization time To in 
Pb+Pb collisions at center of mass energy ^/snn = 2.76 
GeV and impact parameter b = 2.1 fm for the two cases 
of smooth (upper panel) and fluctuating (lower panel) 
initial conditions. We note that the results shown in FigfJJ 
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for the case of fluctuating initial conditions are from one 
of the one hundred different events that are generated for 
this centrality. 



III. THERMAL PROPERTIES OF 
BOTTOMONIA IN QGP 

The potential energy between a pair of heavy quark 
and antiquark is modified in QGP due to the effect of 
color Dcbye screening. The free energy of this heavy 
quark system has been extracted from lattice QCD cal- 
culations, from which the internal energy can then be 
determined from the thermodynamics relation. Whether 
the free energy or the internal energy is more appropri- 
ate for describing the potential energy between a heavy 
quark and antiquark pair in QGP is still controversial 
[301 ] . In this study, we use instead the screened Cornell 
potential (3l| . 



V(r,T) = 



- M (T)r 



« - M (T)r 



(0) 



with a = 0.192 GeV 2 and a = 0.471. The screening 
mass /x(T) depends on temperature, and we use the one 
given in pQCD, i.e., n(T) = (N c /3 + N } /&) 1 / 2 gT, where 
N c and Nf are numbers of colors and light quark flavors, 
respectively, and the coupling constant g is taken to be 
1.87 [22|. Compared to the results from the lattice QCD, 
this potential is close to the free energy around critical 
temperature and becomes more similar to the internal 
energy with increasing temperature. 

Solving the Schrodinger equation with the potential in 
Eq. (|6]) for the bottom quark mass mj, = 4.746 GeV, we 
obtain the dissociation temperatures 681, 285, 190, 257, 
and 185 MeV for the IS, 2S, 3S, IP and 2P bottomonium 
states, respectively [13] • A quarkonium then cannot be 
produced in regions where the temperature is higher than 
its dissociation temperature. 

Even though a quarkonium is produced in less hot re- 
gion, it can be dissociated by scattering with quarks and 
gluons in the QGP. This effect can be quantified by the 
thermal decay width of a bottomonium, 



r(T) =E/ 



d 3 k 



v rel {k)n % {k,T)af s \k,T), (7) 



where i denotes the quarks and gluons in the QGP; rii is 
the number density of parton species i in grand canoni- 
cal ensemble; and v re \ is the relative velocity between the 
scattering bottomonium and parton. For the dissociation 
cross sections of bottomonia af lss , we calculate the m up 
to the next-to-leading order (NLO) in pQCD [H, |33| . 
While in the leading order (LO) a bottomonium is dis- 
sociated by absorbing a thermal gluon, in the NLO it is 
dissociated by the gluon emitted from a quark or gluon 
in the QGP. Fig. [5] shows the thermal decay widths of 
bottomonia as functions of temperature. They are seen 




FIG. 2: (Color online) Thermal decay widths of IS, 2S, IP, 
3S, and 2P state bottomonia shown, respectively, by solid, 
dashed, dotted, dash-dotted, and dash-dot-dotted lines. 



to increase with increasing temperature and diverge at 
their dissociation temperatures. 

Bottomonia are produced with probabilities propor- 
tional to n co \i in Eq. (j3J) for the case of smooth initial 
conditions and Eq. (£[]) for the case of fluctuating initial 
conditions. With its motion isotropically distributed in 
the azimuthal angle 4>, a bottomonium produced at tq 
and moving with velocity v then has the survival proba- 
bility 



S(f ,v) = cxp 



dr r(r, r) 



(8) 



where r = 7% + vr. 

Different from charmonia, bottomonia are hardly pro- 
duced from regeneration in the QGP because of the small 
number of bottom quarks [TtJ ■ Therefore we neglect the 
regeneration effect. 



IV. RESULTS 

Using the thermal properties of bottomonia obtained 
in Sec. lIIIl in the hot dense matter that is described by the 
2+1 ideal hydrodynamic model given in Sec. [HI we have 
calculated their nuclear modification factors in Pb+Pb 
collisions at ^/snn = 2.76 GeV. 

Fig. El shows the Raa of T(1S) as a function of partic- 
ipant number for the smooth (solid line) and fluctuating 
(dashed line) initial conditions together with the experi- 



mental results from the CMS Collaboration [15j. These 
results are obtained with the contributions from Xh(lP), 
Xb(2P), T(2S) and T(35) to T(ls) decays taken to be 
27.1, 10.5, 10.7 and 0.8 %, respectively j34[ and the ve- 
locity of bottomonia determined from the average trans- 
verse momentum of T(IS') (only the transverse momen- 
tum of IS state has been measured.) [l5j]. As discussed 
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FIG. 3: (Color online) Raa of T(1S) as a function of the 
participant number for the smooth (solid line) and fluctu- 
ating (dashed line) initial conditions in Pb+Pb collisions 
at y/SNN — 2.76 GeV. Experimental data are taken from 
Ref. 01. 



FIG. 4: (Color online) Raa of IS (upper lines), 2S (middle 
lines), and 3S (lower lines) bottomonium states as functions 
of transverse momentum for the smooth (solid lines) and fluc- 
tuating (dashed lines) initial conditions in Pb+Pb collisions 
at ^/IWn = 2.76 GeV and b = 2.1 fm. 



in our previous study, most suppression of T(IS') comes 
from the dissociation of its excited states. We note that 
the Raa from both the smooth and fluctuating initial 
conditions are similar to our previous result based on a 
schematic hydrodynamics [l7| , and they are also similar 
to each other except in peripheral collisions. The latter 
is due to the fact that the temperature of the QGP in 
peripheral collisions, which usually does not reach very 
high values in the case of smooth initial conditions, can 
become much higher in the case of fluctuating initial con- 
ditions, leading thus to an enhanced bottomonia dissoci- 
ation. 

Fig. |4] shows the Raa of IS (upper lines), 2S (middle 
lines), and 3S (lower lines) bottomonium states as func- 
tions of transverse momentum in Pb+Pb collisions at 
y?SNN = 2.76 GeV and 6 = 2.1 fm. Solid lines are from 
smooth initial conditions and dashed lines from fluctuat- 
ing initial conditions. The Raa of IP and 2P are similar 
to those of 2S and 3S, respectively. It is seen that the 
Raa of IS state is not changed much by initial fluctua- 
tions as a result of its strong binding and high dissocia- 
tion temperature. The initial fluctuating effect on 2S and 
3S states is, however, not small. Their Raa m the case 
of fluctuating initial conditions are smaller in small pt 
but larger in high px, compared with those in the case of 
smooth initial conditions. As shown in Fig. [TJ nucleon- 
nucleon collisions are more locally concentrated in the 
case of fluctuating initial conditions, resulting in the for- 
mation of hot spots at which there is a relatively larger 
number of binary collisions. Although more bottomonia 
are produced at these hot spots, their survival probabil- 
ity from thermal dissociation decreases unless they have 
enough transverse momentum to escape these regions and 
enhance the so-called leakage effect. As a result, the Raa 



increases more rapidly with transverse momentum in the 
case of fluctuating initial conditions. 
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FIG. 5: (Color online) The ratio of the yield of 2S and 3S bot- 
tomonium states to that of IS state in minimum bias Pb+Pb 
collisions at ^/snn = 2.76 GeV divided by that in p+p colli- 
sions at same energy. Solid and dashed lines are from smooth 
and fluctuating initial conditions, respectively. Experimental 
data are taken from Ref. [lij based on the average transverse 
momentum of IS bottomonium state [T^ |. 

Fig. [5] shows the ratio of the yield of 2S and 3S bot- 
tomonium states to that of IS state in Pb+Pb collisions 
divided by that in p+p collisions at same energy. This 
double ratio has the advantage that the cold nuclear mat- 
ter effect is canceled if they are the same for ground state 
and excited states of bottomonia. The experimental data 
shown in Fig. [5] is 0.31 ±0.03 [l6| at the average transverse 
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momentum of IS state [l5j. The relative yield of 2S and 
3S states in p+p collisions are obtained from Ref. [34 1 . 



Although both smooth and fluctuated initial conditions 
can fit the experimental data, the increase of this dou- 
ble ratio with increasing transverse momentum is steeper 
in the case of fluctuating initial conditions as a result of 
enhanced leakage effect. 



V. SUMMARY 

Bottomonium is a promising particle to probe the 
properties of hot dense matter created in relativistic 
heavy-ion collisions. In this study, we have investigated 
the effect of initial fluctuations in heavy-ion collisions on 
bottomonia suppression. For a more realistic description 
of the expansion dynamics of produced hot matter, a 2+1 
ideal hydrodynamic model was used. The thermal prop- 
erties of bottomonia were obtained from the screened 
Cornell potential and the dissociations of bottomonia by 
thermal partons were calculated up to NLO in pQCD. 
We neglected, however, the small cold nuclear matter ef- 



fect and the regeneration effect. We have found that the 
initial fluctuations hardly affect the survival probability 
of IS state while the effect on excited states is not small. 
It suppresses the survival probability of excited bottomo- 
nia at small transverse momentum but enhances them at 
large momentum, resulting in a survival probability that 
increases more rapidly with transverse momentum than 
in the case of smooth initial conditions. The available 
experimental data on the double ratio of excited states 
to ground state of bottomonium in Pb+Pb collisions to 
p+p collisions at an averaged transverse momentum can, 
however, be described at present with both smooth and 
fluctuating initial conditions. 
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